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Abstract: 

The effect of the environmental interactions on the localization probability 
for a Xe atom trapped in the surface-tip junction of the scanning tunneling 
microscope is studied in the frame of a stochastic, non-linear Liouville equa- 
tion for the density operator. It is shown that the irreversible transfer from 
surface to tip may be explained by thermal decoherence rather than by the 
driving force acting during the application of a voltage pulse. 
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1. Introduction 

Since the first experiments on reversible atom transfer tlie scanning 
tunneling microscope (STM) may be considered as the ideal instrument for 
manipulating atoms or molecules. The bistable operation mode of STM indi- 
cate that in a certain geometry the diffusion barrier on surface is high enough 
to prevent the particle escape from the junction region, and the motion takes 
place along the outer normal to the surface plane in an asymmetric, one- 
dimensional, double-well potential (DWP) p|. However, the mechanism of 
irreversible atom transfer between the potential wells after the application of 
a voltage pulse is not yet completely understood. 

If the initial center of mass (CM) wave function of the atom is an isomeric 
state of the static junction potential, then the barrier crossing appears only 
in special resonance conditions, by quantum coherence oscillations (QCO) 
[0]. The resonances with long oscillation periods require a very fine tun- 
ning of the bias voltage, which cannot be assumed in the present switching 
experiments. Also, during QCO the tunneling has not an exponential law. 
Therefore, is a most interesting issue to study the evolution of the atomic 
CM wave packet dynamics during the voltage pulse, and the decoherence 
effects produced by the coupling with the environment. This coupling may 
damp or completely suppress [§] the QCO. 

The appropriate frame for the treatment of a quantum dynamics deter- 
mined both by the unitary evolution in the Hilbert space and by the change 
in the purity of the states (of the occupation numbers), is provided by the 
Liouville equation for the density operator. However, this equation can be 
solved only for very simple systems, and in a physical situation is necessary 
to find suitable approximations. 

The occupation numbers of the energy levels for a Xe atom in the STM 
junction potential may change in time by the interactions with the electron 
gas. A partial description of this process is given by the rate equations [@], 
but for a realistic treatment is necessary to account also the dynamical effects 
related to the evolution of the non- stationary wave packets, as the QCO. 

In this work, the effect of environmental decoherence on the atom dy- 
namics in the STM junction will be described phenomenologically using a 
modified Liouville equation. This equation is presented in Sec. 2, and is 
applied to calculate the transition rate induced by the thermal noise. In Sec. 
3 are presented numerical results concerning the evolution of the localiza- 
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tion probability for a Xe atom. At zero temperature, the probability density 
during a voltage pulse is determined by the evolution of the quantum wave 
packet in a time-dependent external potential. This evolution is obtained by 
solving numerically the Schrodinger equation, for symmetric triangular and 
trapezoidal pulses with —0.8 V at the peak. The environmental decoherence 
effects are studied using the non-linear Liouville equation for the density op- 
erator, at the temperature T = 4 K and a constant voltage U. Two values 
of U are considered, one non-resonant, U = —0.8 V, and U — 1.141 V, corre- 
sponding to the QCO resonance with the highest frequency. The conclusions 
are summarized in Sec. 4. 



2. Brownian quantum dynamics 



Let us consider a quantum system with the density operator V, inter- 
acting with a classical heat bath of harmonic oscillators via the bilinear 
coupling term ^ 

1=1 

where x is the coupling operator, Cj are constants and the time-dependent 
bath coordinates. The evolution of this mixed classical-quantum system can 
be obtained from the variational equation 

S / dtiY^iQiPi - h,) + Tr[r^^iWtr^ - v^Ho + H,,^,)r]]} = , (2) 
•' 1=1 

where hi = {pf + mjujfqf)/2mi is the classical Hamiltonian for a bath oscil- 
lator, Hq is the Hamiltonian operator for the isolated quantum system, and 
7], are the "square root operators" defined by the Gauss decomposition of 
the density operator, V = rji]'^ p|. In the physical situation of a thermally 
equilibrated bath with infinite heat capacity, Eq. (2) leads to a Brownian 
dynamics of the density operator V, described by the stochastic Liouville 
equation |^ 

ihdtV=[Ho-xm + Mt)),V] . (3) 
Here ^ (t) is a random force with zero mean (the noise) , while 

dTr{Vx) 



,,,)._/;..r,-o^^ (4) 
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is a friction force with the memory function T{t). If << * >> denotes 
the average over the statistical bath ensemble at the temperature T, then 
« ^{t) »= 0, and « ^(t)^(t') »= kBTT{t - t'), by the fluctuation- 
dissipation theorem (FDT). The Brownian evolution of V determined by Eq. 
(3) preserves the "purity" of the initial state, but decoherence may appear 
for the average 

Nt T)r(f\ 

vut)^«v»{t) = Y.^ (5) 

calculated over an ensemble of Nt trajectories T>'^{t) generated with the same 
initial condition, P''(0) = Pq- Of particular interest is the case when the 
initial state of the system is not thermally equilibrated, and Vq =| i^q >< 
ipQ I, with \ ipo > a. pure state. In this case, 

1 Nt 

^^-W = i^ElV^^W><V^^WI (6) 

r=l 

where | ip^it) > is a solution of the modified Schrodinger equation 

ihdt\r>=[Ho-xm + fim\r> ■ (7) 

If the friction may be neglected, then Eq. (3) has the general solution 
with 

P^(t)=re^/o^*'«^(*')^.(*') . (9) 

Here T denotes the time-ordering operator ||^, Ca is the Lie derivative with 
respect to the operator A defined by the commutator, CaB = [A, B], while 

x{t) = e'^^'/'^xe-'""'^" (10) 

is the coupling operator in the interaction representation. 

The evolution determined by the unitary operator exp{—iHQt/h) is the 
same for all the trajectories appearing in Eq. (5), and therefore the ensemble 
average may be written as 

2^a.(t) =e-*^o*/"P,,(t)e^^°*/" . (11) 
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Here T>av{t) denotes the trajectory average of 'D^{t), and can be calculated 
using the FDT after the expansion of the time-ordered exponential in Eq. 
(9). Retaining only the first non- vanishing average, the result is 

T>av{t) = Po - ^ f'dti r dt2r{ti - t2)[x{h), [x{t2),Vo]] . (12) 

n Jo Jo 

This formula can be applied to calculate the rate of the noise-induced tran- 
sitions between the energy eigenstates {| Ek >} of Hq, defined by Hq \ 
Ek >— Ek \ Ek >. If initially the system is in the pure state \ Ei >, then 
T>o —\ Ei >< Ei \, and the rate of the transition | £'j >— >| > is given by 
the asymptotic time-derivative 

A/i = |t-^oo (13) 

of the occupation probability Vf{t) =< Ej \ Vavit) \ Ef >—< Ef \ Vav{t) \ 
Ef >. Using Eq. (12), this probability is 

Vf{t) = 2A;bT^^ /* dti dt2T{h - ^2) cos %(ti - ^2) (14) 
Tl Jo Jo 

where Xfi =< Ef \ x \ Ei > is the matrix element of the coupling operator 
and flfi = {Ef — Ei) /Tl . In the case of Ohmic dissipation with the static fric- 
tion coefficient 7, the memory function is proportional to the delta function, 
T{t) = 27(5(i), and Eq. (13) gives the rate 

2 

A/i = ^ I 1^ ^UbT . (15) 



This was obtained assuming a classical or quasi-classical behavior of the 
environmental degrees of freedom. Therefore, it should provide a good ap- 
proximation when the thermal energy ksT is greater than the transition 
energy, Ti\Vtfi\. 

3. Thermally driven atom tunneling 

The atom dynamics in the STM junction will be treated assuming that 
the CM motion is restricted to the X-axis, normal to the surface, and the 
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potential energy V{x) has an asymmetric, double-well shape. Without exter- 
nal bias, this potential is determined only by the binding interaction energy, 
and the estimate provided in ref. can be well approximated by the fourth 
order polynomial 

Vo{x) = Co + Cix + C2X^ + Csx^ + CiX^ , (16) 

with Co = 0.45 meV, Ci = 0.77 meV/A, C2 = -55.64 meV/A^, C3 = -11.59 
meV/A^, C4 = 44.51 meV/A^. The isomeric minimum of Vo(x) is located 
at Xq = —0.7 A, near surface, the barrier top at Xf, = 0, while the stable 
minimum at Xg = 0.89 A, near the tip. 

At the surface polarization by the voltage U with respect to the tip, the 
dipole interaction energy changes the potential to V{x) = Vo{x) — Ed{x). 
Here E = U/2w is the junction electric field, and d{x) is the dipole moment 
of the Xe atom, 

d{x) = Q.ffx + /io{ Q 3 ^ ^ - Q 3 ^ _ } , (17) 



where Qeff ~ 0.1 e is the average effective charge of Xe HTQ], /^o = 0.3 Debye 
is the induced dipole moment at the surface, w = 2.2 A, and L = 1.56 A. 

If the temperature is zero and there is no dissipation, the atom dynam- 
ics during a voltage pulse U(t) can be obtained by integrating the time- 
dependent Schrodinger equation (TDSE) 

Assuming that U{0) = 0, the initial condition for integration, iIj{x,0) = 
il^o{x), is represented by a Gaussian approximating the isomeric ground state 
of K)(x), 

V;o(x) = (^)i/4e-«>(— 0)^2 (^g) 

TC 



where Co = Muiq/K, M is the Xe mass, and ooo = \J M'^d'^V^/ dx'^\xQ- The 
integration of Eq. (18) was performed numerically in a spatial grid {x^}, 
= 1, A^, by reduction to a Hamilton system of equations. Thus, if Uk{t) = 
Re^ip^Xkyt)) and Vk{t) = Im{%lj{xk,t)) denote the real, respectively the imag- 
inary part of the wave function ip{x,t) at the grid point Xk, then Eq. (18) 
becomes 

d'T'C dT~i 
2Mfc = — 2nvk = -T^ , (20) 



6 



with 

N 



n = Y.Uk{Tu)k + Vk{Tv)k + V{xk){ui + vl) , (21) 



k=l 



(rj. _ , yk+3 + yk-3 ^ yk+2 + yk-2 , ^ Vk+i + Vk-i 49 

2Mdx^^ 90 ^ 20 +^ 2 

The Hamiltonian system of Eq. (20) was defined considering = 321 spatial 
grid points equally spaced by dx = 0.01 A within the interval [xmin,Xmax] = 
[—1.2 A, 2 A]. For a fast integration was used the D02BAF routine of the 
NAG library |ll||, with the time step dt = 6.58 x 10~^ ps. This time step is 



by two orders of magnitude greater than the time step required for the same 
accuracy by the leap-frog method. 

The solution ip{x,t) can be used to calculate the time-evolution of the 
localization probability in the stable well of Vo{x), defined by 

p{t)= / \^{x,t)\'dx . (22) 
Jo 

The results obtained when Qeff = for a triangular voltage pulse of 20 ns are 
presented in Fig. 1(A), and for a trapezoidal pulse of 7 ns, in Fig. 1(C). The 
bias voltage during the pulse is represented in Fig. 1(B) and (D), respectively. 
The localization probability on the tip presented in Fig. 1(A) and 1(C), in- 
creases during the pulse front by sudden jumps at certain time moments 
and t^. These time sequences are different, but the sequences U^, of the 
corresponding bias voltages, defined by = U(t^) and Uj^ = U{t^) with 
U{t) of Fig. 1(B) and (D) are practically the same (e.g. —0.08 V, —0.164 
V, -0.24 V, -0.32 V, -0.4 V, -0.48 V, ...), equally spaced by ~ 0.08 V. 
These jump voltages are also very close to the values of U known to ensure a 
resonant tunneling of the Xe atom between the potential wells by QCO 
Therefore, the jumps are explained by the crossing of the resonances during 
the pulse front, here with arbitrarily chosen slope. 

The features noticed above are present also when Qeff = 0.08 e, as indi- 
cated by the results obtained with a trapezoidal pulse of 7 ns (Fig. 2). The 
increase of p during the pulse front keeps the discontinuous character, but 
the jumps appear at different bias voltages, (e.g. —0.35 V, —0.39 V, —0.43 
V, —0.47 V, .... ), equally spaced by ~ 0.04 V. This result suggests that 
the effective charge represents an additional parameter which may shift the 
values and the spacing of the resonant voltages. In Fig. 2 (A), p becomes 1 
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during the voltage peak, but the residual value after pulse is small. The 
propagation in time up to 20 ns of the non-stationary wave packet created 
by the pulse do not indicate further significant changes in the evolution of p. 

The extension of the upper time limit of the numerical integration to /is 
or ms is not yet possible due to the large amount of computer time required. 
Though, the present results suggest that after the moment when the pulse 
vanishes, the residual probability of localization on the tip is small. 

This behavior was obtained neglecting the dynamical effects produced 
by the coupling between the atom and the surrounding electron gas. The 
average features of the atom-electron interaction are reflected by the jump 
in the junction conductivity at the atom switching, and the potential energy 
term containing the effective charge, assumed above to be fixed. Though, 
the experiments on electromigration in metals show that the effective charge 
is an average quantity which depends on temperature |]T2|. Therefore, it may 
have thermal fluctuations, and for a small value as Qeff ~ 0.1 e, such effects 
could be important. 

In the following, the fluctuating part of Qeff will be treated phenomeno- 
logically, assuming that it can be simulated by an additional dipole interac- 
tion in the Hamiltonian, having a structure close to Hcoup of Eq. (1). 

If the Xe atom adsorbed on the surface is thermally equilibrated, and the 
barrier crossing proceeds by quantum tunneling, then for short times com- 
pared to the recurrence time tr it could be possible to define the average 
transfer rate |1T0 | 

where T* is the effective temperature, Ei are the energies of the isomeric 
levels and Aj the corresponding tunneling rates. However, if the system is 
not equilibrated, and at t = the CM wave function is a pure state ip^, then 
the transfer rate should be defined by 

where Pav{t) =<< p » (t) is the the average localization probability across 
the barrier (x > Xb), in the stable well of V{x), 

PXmax 

Pav{t) = < X \ Vav{t) \ X > dx . (25) 
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According to Eq. (6), the matrix element < x \ Davit) | x > is given by 



< X I Va.{t) \ x>=—J2\<x\ V^^(t) > |2 . (26) 



r=l 



The probabihty amphtude < x \ ip'^{t) >= ip^{x,t) required here can be 
obtained by integrating the equation 

thd.rix, t) = [Ho - xiCit) - 7^^^W{x, t) , (27) 
where < x >r=< ip^\x\ip^ > and 

The average in Eq. (26) was calculated using Nt = 100 solutions ip^{t) of 
Eq. (27) at the environmental temperature T = 4 K and U constant. Each 
solution was obtained considering ^^(t) at the moment t„ = ndt of the form 



Citn) = Rny'^ksT'-) / dt where n = 1,2, 3, ....} is a sequence of Gaussian 
random numbers with mean and variance 1. This choice ensures the discrete 
form of the FDT, « e(^j)e(4) »= 2kBT-f6t^tJdt. 

The static effective charge parameter in d{x) is chosen Qeff = and the 
friction coefficient is 7 = 0.1 /;./ A^. This corresponds to a partial damping 
rate 7/M ~ 5 GHz, small compared to the total estimated rate of vibrational 
relaxation, ~ 30 GHz [|15| . 



At t = the atom is supposed to be in a pure state localized near the 
surface, and the initial condition for integration is the ground state ipo of the 
modified potential 



Vo{x) - ^rf(x) if X < Xb 

Vo{xb) - ^d{xb) + a{x - Xb)2 if X > Xfe 



obtained by replacing the stable well of V{x) by a harmonic term. When U 
decreases near the point when the barrier disappear, this procedure provides a 
better approximation for the isomeric ground state of V{x) than the Gaussian 
wave packet of Eq. (19). 
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The state "^0(2;) was obtained numerically by the Runge-Kutta integration 
of the stationary Schrodinger equation 



[-^^ + Vmod{x)]Mx) = EoMx) , (29) 

where Vmod is defined using a = 0.1. 

The average localization probability on the tip << p » and the average 
energy << Hq »= Tr[HQ << p »] when U = —0.8 V are represented in 
Fig. 3(A) and (B) by solid lines. The time-dependence of << p » can be 
well approximated in this case by the exponential law 

Pfit) = l-e-'^' . (30) 

The rate constant extracted by fit is Ai = 210 MHz, two orders of magnitude 
above the WKB tunneling rate, Xwkb = 6.2 MHz. The WKB rate char- 
acterize the non-dissipative tunneling for times smaller than the recurrence 
time, T/j ~ 2 ps. For longer times, the integration of Eq. (27) with ^ = 0, 
7 = 0, shows that p{t) has the anharmonic oscillatory evolution pictured in 
Fig. 3 (A) by dashed line. These oscillations contain high and low frequency 
components, indicating that the initial wave packet is a superposition of sev- 
eral eigenstates of Hq. During the low frequency oscillation the maximum 
attained by p is pm ~ 0.5 10"'^, while the oscillation amplitude of the average 

position, < X >, is Ar, = {< X >max - < X >min)/2 ~ 10'^ A. 

The low frequency oscillation may be understood assuming that the main 
components of ipo are two eigenstates ipi and ipj of Hq close in energy, but 
localized in the isomeric, respectively in the stable well of V{x). In this case, 
the matrix element \xfi\ is well approximated by ^x/a/p^ = 0.045 A. With 
this value and 7, T used to calculate << p >>, Eq. (15) gives an inelastic 
tunneling rate A/j = 209 MHz, very close to Ai extracted by fit. 

The two-state approximation becomes particularly suited when the poten- 



tial is tuned on a QCO resonance [|13|. At resonance ipo is a linear superposi- 
tion of two quasi-degenerate eigenstates of the Hamiltonian, ip^, ipu each being 
localized with probability 1/2 in either well. Their energies are such 

that the average level spacing 2Tih/TR is much greater than A =| Eu — Ed |, 
and the localization probability of ipo "across the barrier" is 

1 27rt . , 

pQCo{t) = 2 ~ ' (31) 
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oscillating with the period T = 27r/i/A, greater than tr . 

The QCO's of the STM Hamiltonian have an important role in the in- 
crease of the atom transfer probability during the pulse front, but they are 
also very sensitive to the environmental decoherence [jl4[. To study the en- 
vironmental effects on the atomic QCO the bias voltage was fixed at the 
resonant value Ur = —1.141 V when Xe oscillates with the shortest period, 
T = 28.74 ps. Without environment coupling, (7 = 0), the localization 
probability near the tip for the initial wave packet V'o obtained from Eq. 
(29) is pictured in Fig. 4(A) by dashed line. During this oscillation pm ~ 1, 
Ax ^ 0.48 A, and despite the anharmonicities, p{t) is relatively close to 
PQCoif)- At the energy of ipQ the WKB tunneling rate is ~ 50 GHz and 
tr ~ 3 ps. 

For an environmental temperature of T = 4 K, the thermal energy is 
ksT = 0.34 meV, greater than the doublet splitting, A = 0.14 meV, and the 
decoherence effects should be important. The statistical average << p » 
defined by Eq. (25) was calculated using Nt = 100 solutions of Eq. (27), and 
the result is pictured in Fig. 4(A) by solid line. The corresponding average of 
the energy, Tr[ifo << p >>] is shown in Fig. 4(B). 

The average localization probability of Fig.4(A) has damped oscillations 
which can be well approximated by the formula 

pj{t) = 1 - ue-'^' - (1 - u)e-'^' ^^<YT~6^ ' ^^^^ 

The parameters obtained by fit are Ai = 3.6 GHz, A2 = 33 GHz, S = 1.6 ps 
and u = 0.495. 

In a two-state system [tpu, '4'd) coupled to the environment the QCO am- 
plitude is damped with the rate \ud of the intradoublet transitions ipu ^ i'd 
[^], and asymptotically the system arrives in the mixed state of maximum 
entropy 0]. In the present case ipQ is not a superposition of two states only, 
but Eq. (15) with the matrix element \xud\ ~ = 0.48 A gives A^^ = 24 
GHz, close to A2 obtained by fit. 

4. Summary and conclusions 

The atom transfer in the scanning tunneling microscope is a complex 
non-stationary phenomenon, reflecting a dynamical interplay between uni- 
tary evolution in Hilbert space and the environmental decoherence. 
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A phenomenological description of a quantum system interacting bilin- 
early with a classical heat bath of harmonic oscillators is provided by the 
modified Liouville equation presented in Sec. 1. This equation may be 
derived from a variational principle (Eq. (2)), and has a Langevin form, 
containing stochastic and frictional terms. 

The effects of a voltage pulse on the localization probability p for a Xe 
atom prepared initially in a pure state localized on the STM surface was 
investigated by numerically integrating the TDSE (Eq.(18)). In these calcu- 
lations the environmental interactions are neglected, and the voltage pulse is 
assumed of symmetric triangular and trapezoidal shape. The results indicate 
a stepwise increase of p at the moments when the pulse front is near a reso- 
nant bias voltage for the isomeric ground state. The resonant values depend 
on the effective charge, being equally spaced by ~ 0.08 V when Qg// — 0, 
and by ~ 0.04 V when Qeff = 0.08 e. However, the evolution of the non- 
stationary state created by the pulse do not indicate an asymptotic behavior 
characteristic to exponential decay. 

The atom dynamics at finite temperature was investigated in the frame 
of the stochastic, non-hnear Liouville equation (Eq. (3)). The spectrum of 
the environmental noise was assumed to be flat (white noise) , parameterized 
by the static friction coefficient 7. This has a value corresponding to weak 
damping, and the effective temperature was considered the same as the en- 
vironmental temperature, T = 4 K. 

The ensemble average of the energy increases during tunneling, (Fig. 3(B) 
and 4(B)), but within the time interval of 100 ps investigated here the atom 
is not thermalized. The evolution of << p » at resonance {U = —1.141 
V), consists of an incoherent superposition between damped QCO and expo- 
nential decay. In the non-resonant case {U = —0.8 V) the tunneling law is 
close to exponential, with a rate proportional to the product 7T (Eq. (15)). 
This result can provide a basis for understanding the current dependence of 
the atom transfer rate, because the friction coefficient as well as the effective 
temperature should be functions of the electron tunneling current. 
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Figure Captions 



Fig. 1. Atom tunneling at Qeff = without environment coupling, p 
and [/ as a function of time for a triangular pulse of 20 ns (A),(B) and a 
trapezoidal voltage pulse of 7 ns (C),(D). 

Fig. 2. Atom tunneling at Qeff — 0.08 e without environment coupling. 
p (A) and U (B) as a function of time for a trapezoidal voltage pulse of 7 ns. 

Fig. 3. << p » (A, solid) and « Hq » (B) at U = -0.8 V for T = 4 
K, 7 = O.lh/ A^, compared to p at 7 = (A, dash), as a function of time. 

Fig. 4. « p » (A, solid) and « Ho » (B) at U = -1.141 V, for 
T = 4 K, 7 = 0.1^/ A^, compared to p at 7 = (A, dash), as a function of 
time. 
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